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ABSTRACT 

Broad-band (0.8-70 keV) spectra of the persistent X-ray emission from 9 magnetars were obtained 
with Suzaku, including 3 objects in apparent outburst. The soft X-ray component was detected from 
all of them, with a typical blackbody temperature of kT ~ 0.5 keV, while the hard-tail component, 
dominating above ~10 keV, was detected at ~1 mCrab intensity from 7 of them. Therefore, the 
spectrum composed of a soft emission and a hard-tail component may be considered to be a common 
property of magnetars, both in their active and quiescent states. Wide-band spectral analyses revealed 
that the hard-tail component has a 1-60 keV flux, F^, comparable to or even higher than that 
carried by the 1-60 keV soft component, F s . The hardness ratio of these objects, defined as £ = 
Fb/F s , was found to be tightly anti-correlated with their characteristic age r c as £ = (3.3 ± 0.3) x 
(t c /1 ky r )-°- 67 ± - 04 with a correlation coefficient of —0.989, over the range from £ ~ 10 to £ ~ 0.1. 
Magnetars in outburst states were found to lie on the same correlation as relatively quiescent ones. 
This hardness ratio is also positively correlated with their surface magnetic fields with a correlation 
coefficient of 0.873. In addition, the hard-tail component becomes harder towards sources with older 
characteristic ages, with the photon index changing from ~1.7 to ~0.4. 

Subject headings: stars: evolution — stars: magnetars — stars: magnetic field — stars: neutron — 
X-rays: stars 



1. INTRODUCTION 

Soft Gamma Repeaters (SGRs) and Anomalous X-ray 
pulsars (AXPs) are considered to form a distinct class of 
neutron stars, collectively called "magnetars", and are 
thought to have magnetic fields exceeding the critical 
value of -Bqed = 4.4 x 10 13 G (Thompson & Duncan 
1995, 1996; for a recent review see, e.g., Mereghetti 2008; 
Woods & Thompson 2006). Although their radiation, 
mainly emerging in the X-ray frequency, is considered 
to be powered by their huge magnetic energies, little is 
known about how the postulated strong magnetic fields 
are dissipated and converted into the radiation. It is 
hence imperative at this stage to examine their emission 
properties for any empirical relations that provide a more 
unified characterization of these enigmatic objects. 

Through hard X-ray imaging observations with INTE- 
GRAL, some persistently bright magnetars were discov- 
ered to emit not only the well-known soft component 
(with a temperature of kT ~ 0.5 keV), but also a dis- 
tinct hard-tail component which emerges above ~10 keV 
with an extremely hard photon index of Th — 1 (e.g. 
Kuiper et al. 2006; Gotz et al. 2006). The hard-tail com- 
ponent has also been detected from transient magnetars 
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in outburst states, including SGR 0501+4516 (Rea et al. 
2009; Enoto et al. 2010b) and IE 1547.0-5408 (Enoto 
et al. 2010a). This unusual new component, though not 
yet observed from all magnetars, is expected to provide a 
clue to the emission process and the magnetic dissipation 
mechanism in magnetars. 

To perform wide-band spectroscopy of magnetars cov- 
ering these two components, multi-satellite observations 
have generally been employed. However, the hard-tail 
intensity (<10 mCrab) of these objects requires much 
longer exposures than are usually spent in observing their 
soft components. As a result, it has so far been diffi- 
cult to obtain strictly simultaneous wide-band spectra of 
these sources which are more or less time- variable. In 
the present Letter, we report on Suzaku (Mitsuda et al. 
2007) observations of most of the presently known mag- 
netars, studying for the first time soft and hard X-ray 
spectra taken simultaneously. We thus found a strong 
correlation between magnetars' hardness ratio and char- 
acteristic ages. 

2. OBSERVATION AND ANALYSIS 

Since its launch, Suzaku has performed 15 pointed 
observations of 10 magnetars, including 5 targets in a 
Suzaku Key Project in 2009. We analyzed all these 
Suzaku data in a unified way, except for IE 1048.1—5937, 
for which the data are not yet public as of 2009 Novem- 
ber. This included reanalyses of the data from previously 
reported observations; i.e., SGR 1806—20 (Esposito et al. 
2007; Nakagawa et al. 2009), SGR1900+14 (Nakagawa 
et al. 2009), IE 1841-045 (Morii et al. 2008), CXOU 
J164710.2-455216 (Naik et al. 2008), SGR 0501+4516 
(Enoto et al. 2009, 2010b), and IE 1547.0-5408 (Enoto 
et al. 2010a). Table 1 summarizes the observations uti- 
lized in the present Letter. Among them, three sources, 
namely, CXOU J164710. 2-455216, SGR 0501+4516, 
and IE 1547.0—5408, were observed during bursting pe- 
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riods. Our sample, though not complete, covers 9 sources 
out of the ~15 magnetars known to date. 

From each observation, we extracted events for per- 
sistent X-ray emission using the X-ray Imaging Spec- 
trometer (XIS; Koyama et al. 2007) and the Hard X- 
ray Detector (HXD; Takahashi et al. 2007), which cover 
energy ranges of 0.2-12 keV and 10-600 keV, respec- 
tively. The XIS and HXD data were processed with the 
pipeline processing version 2.0 or later, and events were 
discarded if they were acquired in the South Atlantic 
Anomaly or in regions of low cutoff rigidity (<6 GV for 
the XIS and <8 GV for the HXD), or at low Earth eleva- 
tion angles 7 . In addition, a fair number of short bursts, 
detected from SGR 0501+4516, IE 1547.0-5408, and 
SGR 1806—20, were removed by discarding the XIS and 
HXD events which were detected within a time interval 
of a few seconds around each of these bursts (e.g., Enoto 
et al. 2010a). Studies of these burst events are beyond 
the scope of this Letter. 

To create 0.8-10 keV XIS spectra, we accumulated the 
screened data from a source region typically within a 
2'.0 radius of the target centroid, and derived a back- 
ground spectrum from adjacent regions on the same CCD 
chip. Since the background count rates are typically ^5% 
of the source rates, systematic uncertainties of the XIS 
background are negligible. Thanks to the use of appro- 
priate window options, all the observations were free from 
XIS event pile up. 

For each observation, an HXD-PIN spectrum was 
produced in the 10-70 keV range, after subtracting 
non X-ray background (NXB) and Cosmic X-ray back- 
ground (CXB; typically -4% of the NXB). We uti- 
lized NXB events which are simulated by the "tuned" 
model of Fukazawa et al. (2009), and the analytic 
CXB model by Moretti et al. (2009). In the case of 
sources near the Galactic center (i.e., SGR 1806—20 and 
1RXS J170849.0-400910), we further subtracted Galac- 
tic Ridge X-ray Emission (GRXE, Krivonos et al. 2007; 
typically <3% level of the NXB), with a photon index 
fixed at 2.1 (Valinia & Marshall 1998) and normaliza- 
tions specified by observations of nearby blank skies at 
similar Galactic latitudes. Systematic errors in these 
HXD-PIN spectra are dominated by reproducibility of 
the NXB model, which is typically <2% (ler) (Fukazawa 
et al. 2009), while those of the CXB and GRXE are neg- 
ligible. This effect is taken into account by changing the 
PIN-NXB within ±2%. The HXD-PIN field-of-view was 
clear of significant contaminating sources in these obser- 
vations, except for the case of CXOU J164710.2-455216 
(Naik et al. 2008), which was confused with GX 340+0. 
Source confusion was not an issue for the other 8 sources, 
and we hence analyze these 8 sources below. 

3. RESULTS 

Figure 1 shows background-subtracted and pulsc- 
phase-averaged persistent spectra of the 8 magnetars. 
The soft component was clearly observed from all of 
them, while the hard-tail component was successfully de- 
tected from 7 of them at a typical intensity of ~1 mCrab 
(at ^30 keV). Assuming I\ ~ 1, any hard-tail flux from 

7 Detailed criteria are listed in; 
heasarc.nasa.gov/docs/suzaku/proccssing/criteria_hxd.html 
heasarc.nasa.gov/docs/suzaku/processing/critcria_xis.html 



IE 2259+586 was constrained to be less than 8.7 x 10~ 12 
ergs -1 cm~ 2 in the 15-60 keV range, which is consistent 
with the previous result from INTEGRAL (Figure 10 
in Kuiper et al. 2006). These Suzaku spectra recon- 
firm the previous INTEGRAL detections of the hard- 
tail component from some persistent sources, including 
4U 0142+61, 1RXS J170849.0-400910, IE 1841-045, 
SGR 1806-20, and SGR 1900+14 (Mereghetti et al. 
2005; Kuiper et al. 2006; Gotz et al. 2006; den Hartog et 
al. 2007, 2008). The absorbed 2-10 keV and 15-60 keV 
X-ray fluxes are summarized in Table 1. Thus, the hard- 
band flux is comparable to, or sometimes even higher 
than, that in the softer band. 

The unified presentation in Figure 1 yields an in- 
triguing inference that the hard-tail component becomes 
weaker for older objects. To quantify this property, we si- 
multaneously fitted the XIS and the HXD spectra of each 
object with a common spectral model. For the hard-tail 
component, we employed a single power- law for all the 
sources. The soft component of younger magnetars was 
reproduced by a single blackbody (Model A), but older 
sources required some modification to this model due to a 
"soft-tail" excess around 8 keV, which is thought to arise 
via some sort of Comptonzation process (Thompson et 
al. 2002; Lyutikov & Gavriil 2006; Rea et al. 2008). We 
therefore introduced an empirical Comptonized black- 
body model (Model B; Tiengo et al. 2005; Halpern et 
al. 2008; Enoto et al. 2010a), and obtained successful 
fits in most cases. In the case of SGR 0501+4516, nei- 
ther Model A nor B was acceptable, but a successful fit 
was obtained by adding a cooler blackbody to Model B 
(Model C; Enoto et al. 2010b). The XIS spectrum of IE 
1841-045 needed an additional plasma emission compo- 
nent from the surrounding supernova remnant. 

From these fit results, we calculated absorption- 
corrected fluxes for the soft component F s and the hard- 
tail component F^, both defined in the 1-60 keV range. 
These values, given in Table 1, are considered physically 
more meaningful than the raw 2-10 keV and 15-60 keV 
fluxes, because the two components are separated, and 
the absorption is corrected. The lower boundary of 1 
keV was adopted not only for F s but also for F^, because 
sources with young characteristic ages have substantial 
contributions of the hard-tail even in the lowest energy 
range at ~0.8 keV (e.g., SGR 1806-20 in Figure 1). We 
may then define the "hardness ratio" (HR), £ = F^/F s . 
As shown in Table 1, the values of £ range over almost 
two orders of magnitude from 0.2 to 10. 

As visualized in Figure 2, £ exhibits a clear gradient 
on the P-P plane. To assess this implication, we show 
in Figure 3 the values of £ as a function of the charac- 
teristic age t c = P/2Ps and the surface magnetic field 
B s = 3.2 x 10 19 (Ppy/ 2 G, where the errors associated 
with £ include both statistical and systematic errors. The 
latter is dominated by the 2% uncertainties of PIN-NXB 
(§ 2), with additional contributions (by ^1% in £) from 
uncertainties in the XIS vs HXD cross-normalization. 
The figure indicates that £ is indeed negatively (or posi- 
tively) correlated with t c (or B s ). 

After discarding the two upper limits in Figure 3, £ is 
found to be correlated with t c and B s with Pearson's lin- 
ear cross-correlation coefficients of r = —0.989 and 0.873, 
respectively. Thus, £ is very tightly correlated with t c , 
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and to a lesser extent (as evident in Figure 2), with B s . 
In order to evaluate the correlation in a non-parametric 
way, we further performed the Spearman's rank-order 
test, which gives correlation coefficients of r s = —0.963 
and 0.871 for r c and B s , respectively, together with a 
chance probability in both cases of <0.001. Thus, these 
correlations are indeed significant. Fitting through the 
least squares method gave the best-fit linear correlation 
as 

p / \ -0.67±0.04 

-^ = (3.3 + 0.3) x (-£- (1) 



= (0.09 ±0.07) x 



1 kyr 



\Bqed 



1.2±0.2 



(2) 



with Ict uncertainties. Since individual data points have 
rather small errors, we here added constant systematic 
errors to them, so as to make the fits acceptable with 
reduced chi-square values of ~1. Around the £-r c linear 
fit, the data points scatter by no more than a factor of 2 
in £. 

Now that the £-t c and the £,-B s correlations were con- 
firmed to be statistically significant, let us examine them 
for systematic effects. These correlations arc unlikely 
to arise from some selection biases on the P-P plane, 
because these magnetars show a 2-dimensional scatter 
in Figure 2 without any particular correlation between 
P and P. These observations are also free from pos- 
sible uncertainties due to independent variations of F s 
and Fh, since in all sources these two quantities are 
based on strictly simultaneous measurements. In addi- 
tion, multiple observations of the same objects at differ- 
ent epochs gave similar values of £ (i.e., SGR 1806—20, 
SGR 1900+14, and 4U 0142+61 in Figure 3). Therefore, 
the essential features of Figure 2 are not considered to 
depend strongly on the epochs of observations. 

Although the empirical Comptonizcd blackbody model 
(Model B) was employed in calculating F s (and hence 
of 1RXS J170849.0-400910, SGR 0501+45, and 
4U 0142+61, their values of £ change only by a factor of 
1.78, 1.09, and 0.94-1.09, respectively, if their soft com- 
ponent is modeled by two blackbodies that constitute an- 
other popular model. Thus, the implication of Figure 3 
remain unchanged. The £-r c correlation also remains ro- 
bust, even if we change the lower limit of the energy 
band from 1 keV to 0.5 keV. For example, this change 
reduces £ by only 5% and 26% for the cases of SGR 
1806-20 and 4U 0142+61, respectively. Similarly, if we 
raise the upper-limit energy of the hard-tail flux from 60 
keV to 500 keV, the £-t c correlation also persists with 
r = -0.788 and r s = -0.825. Therefore, the £-r c rela- 
tion can be considered to represent a model-independent 
property of our sample of magnetars. Detailed spectral 
analyses of our sample objects will be described else- 
where, together with separate studies of pulsed and im- 
pulsed signals. 

4. DISCUSSION AND CONCLUSION 

Including most of the known SGRs and AXPs, the 
present data set provides the first simultaneous cover- 
age of their soft and hard X-ray spectra (Figure 1). This 
has led to a discovery of the £-r c and £-B s correlations in 
Figure 3. These tight correlations suggest that the values 
of £ in the present objects approximately behave as a one 



parameter family (a single valued function) determined 
uniquely by P and P. Among possible physical quanti- 
ties expressed as simple combinations of P and P, the 
characteristic age r c may be regarded as a likely deter- 
minant of £ (see Figure 2). Then, by regarding r c as an 
indicator of the age in assuming a self-similar spin-down, 
and further noting that our sample covers a significant 
(~50%) fraction of the presently catalogued objects of 
these two classes, the £-r c correlation is considered to 
represent a common spectral evolution of them. Since 
both SGRs and AXPs follow the same correlations, the 
present result provides another support to the increasing 
evidence that these subsets are intrinsically considered 
the same kinds of object, namely "magnetars". 

From the burst-active magnetars SGR 0501+4516 and 
IE 1547.0-5408, we successfully detected the hard-tail 
components. On the £-t c correlation, their £ values gen- 
erally agree with data points of the other magnetars, 
which mainly refer to quiescent or less active states. 
Therefore, the bursting activity of a magnetar is thought 
to cause a significant enhancement in the hard compo- 
nent, roughly to the same intensity as the soft com- 
ponent, which is known to increase by 1-2 orders of 
magnitude (Tarn et al. 2008; Rea et al. 2009; Enoto et 
al. 2010a). The upper limit on the 2nd observation of 
SGR 0501+4516 is consistent with this picture. These 
results further suggest that the basic mechanism of qui- 
escent wide-band emission is retained during increased 
activity. 

In addition to the above interpretation invoking r c , the 
same magnetar £-r c correlation may be interpreted in an 
alternative way. The strong magnetic field B, combined 
with the rotation, would induce an electric field up to 



\E\~mB 



= 2.6 x 10 12 



V 



p 



1/2 



cm 



10 11 ss" 



R 



10 km 



(3) 



where R and are the distance from the center of the 
star and the angular frequency, respectively. Thus, t c 
can be uniquely converted to \E\, as shown at the top of 
Figure 3 (left). Thus, the magnetar £-r c correlation can 

also be regarded as a correlation between £ and \E\. If so, 
the hard-tail component, which should require the pres- 
ence of particles more energetic than ^100 keV, might 
be related to particle acceleration by this induced electric 
field. 

A closer inspection of Figure 1 suggests a hardening of 
the hard-tail component for objects with larger r c . This 
is more clearly shown in Figure 4(a), where is ob- 
served to evolve from ~1.7 to ^0.4. Such a correlation 
was first pointed out in the 1-10 keV band by Marsden 
& White (2001), and recently extended into hard X-rays 
by Kaspi & Boydstun (2010). The latter authors found a 
correlation of with B s or P, with implications similar 
to those of our Figure 4(a). In contrast, they did not find 
any correlations that can be compared with our Figure 3, 
presumably due to lack of simultaneity between the soft 
and hard X-ray measurements. Thanks to the simultane- 
ous observations with Suzaku, we can plot, in Figure4(b), 
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the values of £ as a function of I\. Contrary to normal 
expectations that sources with lower broad-band hard- 
ness ratios should have softer (larger) values of i\, this 
figure reveals an opposite correlation, that older magne- 
tars with weaker hard-tail fluxes (smaller £) in fact have 
harder (smaller) slopes i\. Since this is not affected by 
the selection of energy bands (§3), such a trend should 
imply an intrinsic property of the hard-tail production 
mechanism. 

The present magnetar observations require the hard- 
tail emission mechanism to satisfy both the £-t c corre- 
lation and the hardening of the hard-tail component to- 
wards older objects. The difference of i\ among the 
sources rules out those models which predict a fixed 
common hard-tail slope (e.g., the fast-mode breakdown 
model; Heyl & Hcrnquist 2005). Such a large variation of 
Th with r c is not expected by a thermal bremsstrahlung 
model from the transition layer on the stellar surface, ei- 
ther (Thompson & Beloborodov 2005; Beloborodov & 
Thompson 2007). Similarly, observed Th values have 
not yet been reproduced by existing resonant magnetic 
Compton up-scattering models (Baring & Harding 2007; 
Fernandez & Thompson 2007). It remains to be exam- 
ined whether the results of the present observations can 



Baring, M. G., & Harding, A. K. 2001, ApJ, 547, 929 

Baring, M. G., & Harding, A. K. 2007, Ap&SS, 308, 109 

Beloborodov, A. M., & Thompson, C. 2007, Ap&SS, 308, 631 

Corbel, S., & Eikenberry, S. S. 2004, A&A, 419, 191 

den Hartog, P. R., Kuiper, L., Hermsen, W., Rea, N., Durant, M., 

Stappers, B., Kaspi, V. M., & Dib, R. 2007, Ap&SS, 308, 647 
den Hartog, P. R., Kuiper, L., Hermsen, W., Kaspi, V. M., Dib, 

R., Knodlscder, J., & Gavriil, F. P. 2008, A&A, 489, 245 
Durant, M., & van Kerkwijk, M. H. 2006, ApJ, 650, 1070 
Enoto, T., et al. 2009, ApJ, 693, L122 
Enoto, T., et al. 2010a, PASJ, 62, 475 
Enoto, T., et al. 2010b, ApJ, 715, 665 
Esposito, P., et al. 2007, A&A, 476, 321 
Fernandez, R., & Thompson, C. 2007, ApJ, 660, 615 
Fukazawa, Y., et al. 2009, PASJ, 61, 17 

Gotz, D., Mereghetti, S., Tiengo, A., & Esposito, P. 2006, A&A, 
449, L31 

Halpern, J. P., Gotthclf, E. V., Reynolds, J., Ransom, S. M., & 

Camilo, F. 2008, ApJ, 676, 1178 
Harding, A. K., Baring, M. G, & Gonthier, P. L. 1997, ApJ, 476, 

246 

Harding, A. K., & Lai, D. 2006, Reports on Progress in Physics, 
69, 2631 

Heyl, J. S., & Hernquist, L. 2005, MNRAS, 362, 777 

Hurley, K., Kouveliotou, C, Cline, T., Mazets, E., Golenetskii, S., 

Frederiks, D. D., & van Paradijs, J. 1999, ApJ, 523, L37 
Kaspi, V. M., Gavriil, F. P., Woods, P. M., Jensen, J. B., Roberts, 

M. S. E., & Chakrabarty, D. 2003, ApJ, 588, L93 
Kaspi, V. M., & Boydstun, K. 2010, ApJ, 710, L115 
Koyama, K., et al. 2007, PASJ, 59, 23 

Krivonos, R., Rcvnivtsev, M., Churazov, E., Sazonov, S., 
Grcbenev, S., & Sunyaev, R. 2007, A&A, 463, 957 

Kuiper, L., Hermsen, W., den Hartog, P. R., & Collmar, W. 2006, 
ApJ, 645, 556 

Kulkarni, S. R., Frail, D. A., Kassim, N. E., Murakami, T., & 
Vasisht, G. 1994, Nature, 368, 129 



be explained by any of the existing (or future) theoretical 
explanations (e.g., those based on the "twisted magneto- 
sphere" model; Thompson et al. 2002). 

As an altenative possibility, a photon-splitting pro- 
cess can act to produce the hard-tail emission. At the 
neutron-star surface, sub-MeV photons can be produced 
through electron-positron annihilation or resonant mag- 
netic Compton up-scatterings. Since there is no low 
energy threshold for the photon splitting (Harding & 
Lai 2006), these sub-MeV photons may repeat splitting 
(Harding et al. 1997; Baring & Harding 2001) in the 
strong field, to form the hard continuum downward to- 
wards lower energies. In this case, the higher magnetic 
fields of younger magnetars will allow the photons to 
repeat splitting toward lower energies, thus making I\ 
larger in agreement with our figure. 
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TABLE 1 

Log of Suzaku observations of magnetars. 



Name" P b P 

10 -11 10 
(s) (..-'KG) 



11 nn 14 



Obs. ID 



(kyr) 



Date 



yymmdd 



Exp." 



Model 



(kg) (xVd.o.f) 



Abs. F x s 

P15-60/ P2-10 
15-60/2-10 



Unabs. F x n 

Ph/P B 
1-60/1-60 



HR' 

~-F h /F s 



rv 



1806-20 7.556 54.9 21 0.22 

1900+14 5.169 7.78 6.4 1.1 
0501+45 5.762 0.5 1.7 18 



Soft Gamma-ray Repeater 

401092010 06-09-09 55.4 A(468. 9/434) 

401021010 07-03-30 16.5 A(185. 1/182) 

402094010 07-10-14 50.1 A(239. 3/240) 

401022010 06-04-01 14.4 A(145. 5/156) 

404077010 09-04-26 42.1 A(433. 6/393) 

903002010 08-08-26 33.5 0(151.2/149) 

404078010 09-08-17 43.0 B(191. 7/154) 
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Anomalous X-ray Pulsars 

903006010 09-01-28 33.5 A(298. 9/278) 105(7)/57(1) 

401100010 06-04-19 63.5 B(462. 9/262) 45(4)/26(2) 

404080010 09-08-23 60.9 B(94.6/98) 27(6)/28.8(6) 

402013010 07-08-13 101.6 B(434. 4/408) 36(3)/64.4(2) 

404079010 09-08-12 107.4 B(310.6/301) 30(4)/61.4(3) 

901002010 06-09-23 35.1 B(333.2/294) — /23.4(3) 

404076010 09-05-25 103.4 B(222. 1/213) <8. 7/15.8(1) 
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Abbreviated names represent SGR 1900+14, SGR 1806-20, SGR 0501+4516, IE 1547.0-5408, IE 2259+586, 4U 0142+61, 
CXO J164710. 2-455216, 1RXS J170849.0-400910, and IE 1841-045. 
b P and P values are from the McGill SGR/AXP Online Catalog (http://www. physics. mcgill.ca/~pulsar/magnetar/main. html). 
c Surfacc magnetic field strength estimated as B s = 3.2 X 10 19 (PP) 1//2 G. 
d Characteristic ages estimated as r c = P/2Ps. 
° Effective Exposures of the HXD. 

f Model A: a blackbody plus a power-law model (BB+PL); Model B: a Comptonized blackbody plus the power-law (CBB+PL); Model C: 
a blackbody and the Comptonized blackbody plus the power-law (BB+CBB+PL). 

g Absorbed X-ray fluxes (10 — 12 ergs -1 cm~ 2 ) in the 15-60 keV and 2—10 keV bands. Quoted errors arc statistical-only at the 90% 
confidence level. 

h Unabsorbcd X-ray fluxes (10 -12 ergs -1 cm -2 ) of the 1—60 keV soft component and the 1-60 keV hard-tail component. Statistical errors 
(lcr) arc included. Systematic errors (lcr) are also included in the hard-tail component, calculated from the uncertainties in the PIN-NXB 
and in the cross-normalization between the XIS and HXD-PIN. 

'Hardness ratio £, calculated from the two unabsorbed fluxes. Quoted errors are statistical and systematic lcr level. 

J Photon index of the hard-tail component. Quoted errors are statistical and systematic lcr level. If the upper and lower errors are different, 
average errors are shown. 
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Energy (keV) 

Fig. 1. — Background-subtracted vF v spectra of the persistent emission of the magnetars, shown after eliminating the instrumental 
responses. Interstellar absorption is included. Individual spectra are shown with offsets, and are arranged in order of increasing characteristic 
age from top (young) to bottom (old). Blue horizontal lines indicate a 1 mCrab intensity. Green, red, and blue lines represents the soft 
component, the hard-tail component, and the SNR contamination (including line emission) in the IE 1841-045 spectrum, respectively. If a 
source was observed more than once, one observation is shown. The GSO data of IE 1547.0—5408 are included after Enoto et al. (2010a). 
The full names of the sources, which are shown here as abbreviations, are listed in a footnote to Table 1. 



Broad-band study with Suzaku of the magnetar class 
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Fig. 2. — A portion of the P-P diagram of pulsars (Manchester et al. 2005; http://www.atnf.csiro.au/research/pulsar/psrcat/). Sizes 
of the red circles correspond to the £ values in the present sample. Blue and green lines represent constant values of the magnetic field 
strength and the characteristic age, respectively. Small grey circles represent magnetars with upper limits on the hard-tail component, or 
those not observed with Suzaku. Ordinary radio pulsars are indicated by crosses. 
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Fig. 3. — (left) A correlation between the HR £ and the characteristic age r c . Quoted errors are statistical plus systematic 1<t (Table 1). 
Upper limits are also shown for SGR 0501+4516 (2nd observation) and IE 2259+586. Green solid and dashed lines represent the best fit of 
equation (1) and their boundaries shifted by a factor of two, respectively. SGRs and AXPs are shown in red and blue, respectively, while 
burst-active sources are shown with star symbols. (Right) The same as the left panel, but plotted as a function of the surface magnetic 
field. 
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Fig. 4. — (a) Photon indices Tj, of the hard-tail component of the present sample, plotted as a function of the characteristic age. Quoted 
errors are statistical and systematic lcr. Colors and symbols are the same as in Figure 3. Green squares and triangles represent photon 
indexes of pulsed and total fluxes, respectively, after Kaspi & Boydstun (2010). (b) Relation between and f. 



